arXiv:1507.08002vl [astro-ph.SR] 29Jul2015 


MNRAS 000.[TKTT|(2015) 


Preprint 30 July 2015 


Compiled using MNRAS IATpX style file v3.0 


New members of the TW Hydrae Association and two accreting 
M-dwarfs in Scorpius-Centaurus 


Simon J. Murphy 1,2 *, Warrick A. Lawson 2 and Joao Bento 1 

1 Research School of Astronomy & Astrophysics, Australian National University, Canberra, ACT 2611, Australia 

2 School of Physical, Environmental and Mathematical Sciences, University of New South Wales Canberra, ACT 2600, Australia 


30 July 2015 


ABSTRACT 

We report the serendipitous discovery of several young mid-M stars found during a search for 
new members of the 30-40 Myr-old Octans Association. Only one of the stars may be consid¬ 
ered a possible Octans(-Near) member. However, two stars have proper motions, kinematic 
distances, radial velocities, photometry and Li I A6708 measurements consistent with mem¬ 
bership in the 8-10 Myr-old TW Hydrae Association. Another may be an outlying member 
of TW Hydrae but has a velocity similar to that predicted by membership in Octans. We also 
identify two new lithium-rich members of the neighbouring Scorpius-Centaurus OB Associ¬ 
ation (Sco-Cen). Both exhibit large 12 and 22 |xm excesses and strong, variable Ha emission 
which we attribute to accretion from circumstellar discs. Such stars are thought to be incred¬ 
ibly rare at the ~16 Myr median age of Sco-Cen and they join only one other confirmed M- 
type and three higher-mass accretors outside of Upper Scorpius. The serendipitous discovery 
of two accreting stars hosting large quantities of circumstellar material may be indicative of 
a sizeable age spread in Sco-Cen, or further evidence that disc dispersal and planet formation 
time-scales are longer around lower-mass stars. To aid future studies of Sco-Cen we also pro¬ 
vide a newly-compiled catalogue of 305 early-type Hipparcos members with spectroscopic 
radial velocities sourced from the literature. 

Key words: open clusters and associations: individual: TW Hydrae, Scorpius-Centaurus (Sco 
OB2) - stars: pre-main sequence - stars: formation - stars: low-mass 


1 INTRODUCTION 

Nearby stars with ages of 5-20 Myr are ideal laboratories in which 
to investigate the later stages of star and planet formation. It is on 
these time-scales that the gas-rich discs which feed giant planet 
formation and circumstellar accretion dissipate and are replaced 
by dusty debris discs replenished by the collisions of rocky bodies 
( |Wyatt|2008| |Williams & Cieza|20fl| . Such stars are also prime 
targets for direct imaging observations of young, bright exoplanets 
and substellar companions (e.g. Marois et al. 2008; Lagran ge et al.| 
|2010| >. The discovery of young moving groups in the solar neigh¬ 
bourhood over the last two decades (see reviews by |Zuckerman &| 
|Song|2004||Torres et al.|2008) has enabled the study of circumstel¬ 
lar material and processes in this important age range at high spatial 
resolution and sensitivity across a wide range of stellar masses. 

Foremost amongst the nearby moving groups is the retinue 
of coeval pre-main sequence stars (age 8-10 Myr; Weinberger, | 


|Zuckerman et al.|2001||Song^Zuckerman&Bessell 2003M>ut also 


|Anglada-Escude & Boss||20T3j |Ducourant et al.||20 1 4^ within a 
few tens of degrees and comoving with the T Tauri star TW Hydra 


(jRucinski & Krautter 1983 

. The so-called TW Hydrae Associa- 

tion (TWA; de la Reza et al. 

1989 Kastneretal. 1997, Torres et al. 

2008; Schneider, Melis & Song 2012a and references therein) now 


includes approximately 30 systems, spanning spectral types of A0 
through to brown dwarfs. Its members were identified by a vari¬ 
ety of means, most commonly as common proper motion sources 
detected in X-rays by the ROSAT mission (e.g. [Webb et al.|1999[ 

* E-mail: simon.murphy@anu.edu.au 


through near-IR photometry JGizisJ2002J Looperetah|2007 


circumstellar disc excesses (|Gregorio-Hetem et al.|1992[ Schneider| 


2010 ). 


|et al,|2 012bl and activity ( [Shkolnik et al,|20 1 1 1 . The current mem¬ 
bership of TWA is dominated by K and M-type stars, although there 
appears to be a dearth of members later than ~M3 ( [Shkolnik et al.| 
[2011[ ). This is partly a consequence of the heavy reliance of pre¬ 
vious searches on the ROSAT All Sky Survey 1 Voge s~et al.|1999} >, 
which, for the saturated coronal X-ray emission commonly seen in 
young stars (Lx /Lbol ~ — 3), had a limiting spectral type of around 
M2 at the ~60 pc mean distance of TWA. Recent efforts exploit¬ 
ing deeper UV photometry from the GALEX mission l | Mart in et al.| 
|2005p have proven successful at identifying young stars from their 
enhanced chromospheric activity [Shkolnik et al.|20lT[|Rodriguez| 
|et al.||20TT] |Murphy & Lawson 2015| ), especially in conjunction 

with modern photometric and proper motion surveys. 

Immediately south of TWA lies the Scorpius-Centaurus OB 
Association (Sco-Cen, Sco OB2; Blaauw 1964; de Zeeuw et al. 
|1999| |Preibisch & Mamajek| |2008j >, the closest site of recent 
massive star formation and the dominant population of pre-main 
sequence stars in the solar neighbourhood. The association has 
been implicated in the formation of several nearby young moving 
groups, including TWA (Mamajek & Feigelson 2001 ; [Femandez,| 

|Figueras & T orra 2008 : , |Murphy, Lawson & Bessell|2013| . Sco-Cen 

has traditionally been divided into three subgroups (Fig. |TJ, each 
with subtly different mean distances, ages and velocities: Upper 
Scorpius (US, 145 pc), Upper Centaurus Lupus (UCL, 140 pc) and 
Lower Centaurus Crux (LCC 120 pc; |de Zeeuw et al.|1999j l. The 
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Figure 1. Galactic positions and proper motions of the eight lithium-rich stars identified in this work (black points). Also plotted are proposed Sco-Cen (blue 
dots; |de Zeeuw et a1. 1 1 999| and TWA members (red squares; Ducourant et al. 2014). Open TWA symbols denote non-members, most of which are probably 
background LCC stars. The star near l = 280 deg with the large proper motion is TWA 22, a likely j3 Pictoris member (Teixeira et al. 20091. For reference, the 
well-known Southern Cross asterism is also plotted at / = 300 deg. 


subgroups have median ages of approximately 10 Myr, 16 Myr and 
17 Myr, respectively ((Mamajek, Meyer & Liebert||2002[ [Pecauh] 
|Mamajek & Bubar||2012}. While the full mass function of US is 


well known (e.g. Preibisch^£t^alJlj98 yLodieu, Dobbie & Harm] 

|bly||201l| |Rizzuto, Ireland & Kraus 2015 i, the older and more 
southern UCL and LCC subgroups are under-studied below 1 Mq. 
Only ~100 K and M-type members are known, mostly detected by 
ROSAT i Mamajek et al.|2002[|Preibisch & Mamajek|2008[[Song, 


Zuckerman & Bessell 2012) but also more recently by GALEX ( Ro 
driguez‘etaT|20TT|>. Any reasonable mass function predicts several 


thousand low-mass members awaiting discovery across Sco-Cen. 

In this contribution we describe the serendipitous discovery 
and spectroscopic confirmation of one existing and up to three 
new M-type members of TWA, as well as two previously unknown 
members of Sco-Cen exhibiting signs of accretion from circumstel- 
lar discs. These six stars, and another two spectroscopically-young 
stars not yet attributable to any known group, were identified dur¬ 
ing an ongoing search for members of the 30-40 Myr-old Octans 
Association ( jTorres et al.[20()8~] |Murphy & Lawson|2015} . In the 
following section we briefly describe the candidate selection pro¬ 
cedures and observations, deferring a full description of the Octans 
survey to an upcoming work (Murphy et al. 2015, in preparation). 


2 CANDIDATE SELECTION AND OBSERVATIONS 

The aim of our photometric and proper motion survey was to 
identify new Octans members on the cool side of its lithium de¬ 
pletion boundary (LDB) in order to better refine the age of the 
association. Candidates for spectroscopy were selected from the 
SPM4 proper motion catalogue i Girard et al.||20fT| l in a sim¬ 
ilar manner to jMurphy & Lawson {2015 1 


with initial cuts of 
Uspm4 < 20, V — K > 5, J — H < 0.7 and proper motion errors of 
less than 10 masyr -1 . From this list we then selected stars whose 
proper motions agreed with th e Octans space motion, (U ,V,W) = 
(—13.7,—4.8, —10.9) kms -1 (Murphy & Lawson 2015), at better 
than 2oy,, where Oj; is the total proper motion uncertainty. Using 
the resulting kinematic distances we also required each candidate 
lie within the spatial bounds of the Octans/Octans-Near ( |Zucker-| 
|man et al.|2013} complex, —100 < X < 150 pc, —150 < Y < 50 pc 
and —100 < Z < 50 pc. To further limit the sample, we finally 


required a match against near-UV photometry from the GALEX 
mission. Following these cuts, 369 candidates remained over the 
SPM4 footprint (5 < —20 deg). After visual inspection to eliminate 
galaxies and objects with obviously erroneous photometry, spectro¬ 
scopic observations concentrated on the reddest stars lying above a 
Pleiades isochrone ( Stauffer et al. 2007 1 with the best proper mo¬ 
tion matches and near-UV excesses {Rodriguez et al.|20lift . 

We obtained spectroscopic observations of 25 high-priority 
candidates and a selection of velocity and spectral type standards 
using the Wide Field Spectrograph (WiFeS; |Dopita et al.|2007| on 
the Australian National University 2.3-m telescope at Siding Spring 
Observatory during 2015 March 6-11. The instrument set up, data 
reduction and analysis, including the derivation of spectral types 
and radial velocities, was the same as that described in |Murphy &| 
|Lawson|{2015| > and we refer the reader to that work for more infor¬ 
mation. All spectra were flux calibrated using observations of the 
white dwarf EG 21 {Bessell|1999| l. Further WiFeS observations of 
two interesting candidates (2M1239-5702 and 2M1422-3623, see 
Section |3~2) were taken with the same instrument settings between 
2015 March 30 and April 6 as part of another programing] 


3 RESULTS 

Only 4/25 candidates had radial velocities within 5 km s _1 of those 
predicted by Octans membership. The relationship of these stars 
(only one of which is lithium-rich) to Octans will be discussed in a 
future work. Eight candidates showed strong Li I A 6708 absorption 
(EW(Li) > 100 mA). These stars are listed in Table [I] and their 
WiFeS R 7000 spectra are shown in Fig. [2] The detection of lithium 
is an indicator of youth in M-type stars, but with a strong stellar 
mass dependence {Soderblom|2010) . At ages of 10-20 Myr lithium 
absorption disappears at early and mid-M spectral types only to 
reappear at lower masses as one crosses the LDB (e.g. |Basri, Marcy| 
& Graham 1996). With spectral types of ~M5, the eight stars in 


Table 


are demonstrably younger than the Pleiades, (LDB M6.5; 


[Stauffer, Schultz & Kirkpatrick) 1998j ), and are likely no older than 

1 Reduced spectra of the lithium-rich stars and machine-readable versions 
of Table Al can be found on the FIGSHARE service at the following DOI: 
http: //dx. doi . org/10.6084/m9 . f igshare . 1491506 
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Table 1. Lithium-rich candidates observed with 2.3-m/WiFeS during 2015 March-April. Spectral types are estimated from V — K s colours and comparison to 
standards. V magnitudes are from SPM4, except for 2M1058 -2346 and 2M1202—3328 which are from APASS-DR6. Photometry for 2M1028-2830,2M1200- 
3405 and 2M1422-3623 are derived from photographic plates and will be of poorer quality than SPM4 CCD photometry. 2M1200-3405 has a broadened cross 
correlation function (CCF) which may indicate binarity (see text). Candidates displaying emission lines in their spectra are noted (see Table[3j. 


2MASS designation 
(72000) 

Spectral Type 
(±0.5) 

V 

(mag) 

K, 

(mag) 

EW(Li) 
(±50 mA) 

EW(Ha) 
(±1 A) 

vio(Ha) 
(±10 kms~*) 

RV 

(±2 kms -1 ) 

Notes 

2MASS J10055826-2131142 

M4.9 

16.72 

10.70 

630 

-11.5 

172 

11.0 

- 

2MASS J10284580-2830374 

M4.9 

16.03 

10.03 

600 

-10.9 

162 

20.0 

TWA 34; [OI] 

2MASS J10491880-2509235 

M4.9 

17.08 

11.08 

800 

-11.5 

186 

19.0 

- 

2MASS J10585054—2346206 

M5.1 

15.58 

9.43 

700 

-11.3 

151 

2.0 

- 

2MASS J12002750-3405371 

M4.7 

14.55 

8.72 

650 

-9.9 

312 

11.0 

Broad CCF 

2MASS J12023799-3328402 

M4.8 

15.77 

9.85 

700 

-9.5 

169 

6.0 

- 

2MASS J12392312-5702400 

M5.0 

17.35 

11.25 

650 

[-63, -27] 

[238,331] 

16.0 

Hef 

2MASS J14224891-3623009 

M5.0 

17.31 

11.27 

600 

[-91, -33] 

[236, 341] 

9.0 

He I, [OI], Na f 




Figure 2. 2.3-m/WiFeS R 7000 spectra of the lithium-rich stars identified in this work. The full 5300-7050 A spectra are plotted in the left panel and the 
shaded region around Ha and Li I A 6708 is enlarged on the right. All spectra are normalised over 6100-6150 A. The broader spectral features of 2M1200- 
3405 (especially around Li I) are readily apparent. Only first epoch spectra of 2M1422-3623 and 2M1239-5702 are plotted (see Table[3j. 


the nearby ~20 Myr /3 Pictoris Association (LDB M4.5; Binks & 
|jeffries|2014j l. The sky positions and proper motions of the lithium- 
rich stars are plotted in Fig.[I]in Galactic coordinates. All either lie 
within or near TWA and the adjacent Sco-Cen. 


3.1 New members of the TW Hydrae association 

2M1058-2346, 2M1200-3405 and 2M1202-3328 have positions 
and proper motions similar to nearby TWA members (Fig. [TJ. The 
latter two stars lie within 2 deg of the confirmed Ml member 
TWA 23 jSong et al.|2 003|l and all three have strong lithium ab¬ 
sorption similar to other M-type members. Shk olnik et al.| |20TT| 
identified 2M1200-3405 as a UV-bright TWA candidate but did not 


report any Ha emission in their Magellan/MagE spectrum. Based 
on this non-detection they disregarded 2M1200-3405 as a possible 
member and did not obtain a higher-resolution spectrum or radial 
velocity. However, subsequent examination of the raw MagE frame 
shows clear Ha emission and that the star was inadvertently ex¬ 
cluded as a likely member (E. Shkolnik, personal communication). 
The WiFeS spectrum of 2M1200-3405 plotted in Fig.[2]has strong 
Ha emission (EW = —9.9 A) typical of an active mid-M dwarf. 
The broad velocity width (vio = 312 kms -1 ) may be the result of 
spectroscopic binarity (see Section [3. 1.1 We find 2M1200-3405 
has a GALEX NUV — J excess consistent with other nearby young 
stars l |Rodriguez et al.|201l| l, while 2M1058-2346 and 2M1202- 
3328 are borderline UV-excess sources. 
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To support the membership of 2M1058-2346, 2M1200-3405 
and 2M1202 -3328 in TWA we use the BANYAN II web service 
(version 1.3; Malo et al7||2013 Gagne et al.||2014p BANYAN is 
a Bayesian analysis tool which compares a star’s Galactic posi¬ 
tion ( XYZ ) and space velocity ( UVW) to the observed distribu¬ 
tions of young groups and the field. No photometric information is 
used. It utilises a naive Bayesian classifier to calculate membership 
probabilities, marginalising over the unknown distance and taking 
into account the expected populations (prior probabilities) of each 
group and the field, as well as uncertainties on the observables. 
Within BANYAN we adopt the ‘young’ field hypothesis (age <1 
Gyr) given the clear youth of the lithium-rich stars. Using only their 
sky positions and SPM4 proper motions, BANYAN returns TWA 
membership probabilities of 83, 99 and 98 per cent for 2M1058- 
2346, 2M1200-3405 and 2M1202-3328, respectively, with the bal¬ 
ance of membership probabilities going to the field hypothesis. The 
predicted radial velocities of 10.3 ± 1.4 and 9.9 ± 1.5 kms -1 for 
2M1200-3405 and 2M1202-3328 agree well with the WiFeS val¬ 
ues in Table[l] further supporting the stars’ membership in TW40 
Including these observations in BANYAN gives membership prob¬ 
abilities of 99 and 97 per cent, with implied distances of ~60 pc 
for both stars. The resultant space motions of 2M1200-3405 and 
2M1202-3328 are listed in Tableland plotted in Fig.[3]with con¬ 
firmed TWA members. Their velocities are in excellent agreement 
with the other members, as expected. The mean space motion of the 
24 members in this diagram is ( U,V,W ) = (-10.8,-18.3 — 5.0)3= 
(2.0,1.7,1.6) kms" 1 (let variation). BANYAN uses a different set 
of 18 members in its TWA model, whose principle spatial and ve¬ 
locity axes are not aligned to standard Galactic coordinates (see 
|Gagne et al.|2014^ However, both samples have almost identical 
mean velocities (Av < 1 km s" 1 ). 

On the contrary, the WiFeS radial velocity of 2M105 8-2346 
(2 =t 2 kms -1 ) disagrees by 4 a with that predicted by BANYAN 
for membership in TWA. However, it is the only lithium-rich star 
to have a space motion consistent with membership in Octans 
(Fig. i RV 0ct = —2.2 kms ') and the only one associated with 
a ROSAT source (1RXS J105849.9-234623; [Voges et al.[[T999 j. 
We estimate Fx ~ 2 x 10~ 13 erg s" 1 cm" 2 3 (near the ROSAT de¬ 
tection limit) and logL^/Lbol ~ — 3.2 using the relations in 


Flem- 


|mg, Schmitt & Giampapa| ([l995) and Pecaut & Mamajek (2013| >. 
Detection by ROSAT at this spectral type immediately implies the 
star is closer than the majority of TWA members. BANYAN pre¬ 
dicts a distance of 41 ±4 pc, similar to that implied by member¬ 
ship in Octans (38 pc), placing it near the periphery of the ‘Octans- 
Near’ group ( Zuckerman et al.j2013||Murphy & Lawsonf2015} . If 
Octans-Near is the same age as Octans and similar to the Tucana- 
Horologium Association (30-40 Myr; see discussion in |Murphy &| 
|Lawson|2015) , the strong lithium absorption observed in 2M1058- 
2346 (EW = 700 mA) implies it is on the cool side of the Octans 
LDB, the first such member to be discovered (see Fig. [4}. Alterna¬ 
tively, given its sky position and proper motion vector it is possible 
2M105 8-2346 is a spectroscopic binary member of TWA or that 
its single-epoch WiFeS radial velocity is erroneous. In this case its 
lower membership probability in BANYAN could be due to its lo¬ 
cation near the edge of the TWA spatial ellipsoid. Further radial 
velocity measurements are necessary to distinguish between these 
two scenarios. 


2 http://www.astro.umontreal.ca/~gagne/banyanll.php 

3 2M1202-3328 was observed with the ESO 2.2-m/FEROS (programme 
ID 090.C-0200A) on 2013 February 17 as part of the GALEX Nearby 
Young-Star Survey (GALNYSS; |Rodriguez et al .] |2011| |20 ] 3} . The 
pipeline-reduced spectrum is heavily contaminated by residual sky emission 
and strong noise spikes. By combining line velocities for Ha, H/3 and the 
Na I 78200 absorption doublet we derive a mean radial velocity of 8.3 ± 1.7 
(s.d.) kms" 1 , in good agreement with the WiFeS and BANYAN values. 



U (km s-') 

Figure 3. Candidate space motions (see Table [2), coloured by their pro¬ 
posed memberships: 2M1200-3405 and 2M1202-3328 (TWA, red points), 
2M 1028-2830/TWA 34 (TWA, purple star), 2M105 8-2346 (Octans/TWA?, 
green diamond), 2M1239-5702 and 2M1422-3623 (Sco-Cen, orange trian¬ 
gles) and 2M1005-2131 and 2M1049-2509 (unknown, blue squares). Also 
plotted are confirmed TWA members with trigonometric parallaxes (black 
points, with ±2 ct range shaded;|&ucourant et al.|2QI4f. 
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Figure 4. Li I A6708 EW of 2M1058-2346 compared to members of TWA 
ISchneid er et al.[2012a|, T uc-Hor |da Silva et al.|2009||Kraus et al.|2014) 

and Octans 7Mmphy^^awson|2Q15t . The 30 40 MyrA)ld Tuc-Hor lithium 

depletion boundary (LDB) is clearly visible at V — K s w 6 (dashed line). 


As a final check, we plot in Fig. [5] the My versus V — K s 
colour magnitude diagram (CMD) for the lithium-rich stars, to¬ 
gether with empirical sequences for TWA, /3 Pictoris and the 
Pleiades. 2M1058-2346 and 2M1202-3328 have CCD-derived 
APASS ( jHenden et al.|20f2| > V magnitudes which we adopt over 
the SPM4 values (but which agree within 0.2 mag), while for 
2M1028-2830, 2M1200-3405 and 2M1422-3623 the SPM4 V 
magnitude was measured from photographic plates and should 
be treated with caution. No APASS data exist for these stars. 
Both 2M1200-3405 and 2M1202-3328 fall close to the mean lo¬ 
cus of TWA members l|Riedel et al.|j20lT|. In conjunction with 
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Table 2. Kinematics and memberships of the lithium-rich stars. All proper motions are from SPM4. Distances are from BANYAN (BAN) or purely kinematic 
estimates from the projected group velocity (kin). The latter are assumed to have a 20 per cent parallax error (10 per cent for 2M1239—5702). UVW velocities 
define a right-handed triad with U positive towards the Galactic centre. Uncertainties were calculated following the prescription of |Johnson & Soderblom| 
(1987) using the uncertainties in proper motion, radial velocity and parallax. The space motion of 2M1028-2830 (also see Fig. [3) assumes the higher precision 
radial velocity (13.3 ± 0.1 kins 1 ) obtained by Rodriguez et al. from ALMA observations (see text). 


2MASS designation 

Ha COS 8 

bs 

Dist. 

Ref. 

u 

V 

W 

Membership 

(72000) 

(mas yr -1 ) 

(mas yr -1 ) 

(pc) 


(kms~*) 

(kms~*) 

(km s -1 ) 


2MASS J10055826-2131142 

—41.7±4.6 

—13.4±4.2 

75+ 19 
' -12 

kin. 

—10.8±2.4 

— 14.1 ±2.1 

—6.9±2.9 

? 

2MASS J10284580-2830374 

—68.6±2.7 

— 11.4 ± 2.5 

47 ±6 

BAN. 

—11.5± 1.5 

-16.3±0.6 

-4.3 ±1.4 

TWA (34) 

2MASS J10491880-2509235 

—39.0±2.0 

—7.9 ±1.8 

95±f* 

kin. 

—12.8±2.7 

- 22.8 ± 2.2 

— 1.1 ±2.4 

? 

2MASS J10585054—2346206 

—98.3 ± 1.9 

—20.6± 1.7 

381 \ 

kin. 

—13.2±2.7 

— 8 . 6 ± 2.2 

-9.1 ±2.3 

Oct-Near/TWA? 

2MASS J12002750-3405371 

—55.6±3.2 

-19.2±3.0 

61 ±7 

BAN. 

—8.5 ±1.7 

—18.2±2.0 

—2.9± 1.5 

TWA (35) 

2MASS J12023799-3328402 

—58.4±3.8 

—17.3 ±3.5 

62 ±7 

BAN. 

—11.2 ± 1.9 

—14.4±2.0 

—4.8± 1.6 

TWA (36) 

2MASS J12392312-5702400 

—32.0±3.5 

—10.4±3.7 

11911? 

kin. 

-6.6±2.5 

—23.4±2.2 

-5.1 ±2.2 

Sco-Cen (LCC) 

2MASS J14224891-3623009 

—17.0±6.3 

—25.6±5.8 

142lg 

kin. 

-0.1 ±3.5 

—21.4±4.8 

—7.2±4.3 

Sco-Cen (UCL) 


MO M2 M4 M5 M6 



Figure 5. Candidate colour-magnitude diagram. The grey point assumes 
2M1200-3405 is a corrected equal-brightness system (see text). Also plot¬ 
ted are members of TWA (red squares; Ducourant et al. 20141 and Tuc-Hor 
(blue dots;|Kraus et al.|2014). Lines show fits to TWA (red) and j 8 Pic (or¬ 
ange) members from|Riedel et al.||201 1| and the Pleiades sequence (Stauf-| 
|fer et al.|2007) . The maximal error bar represents a 20 per cent parallax 
error and 0.1 mag photometric errors. 


their kinematic and spatial match from BANYAN, the CMD place¬ 
ment, strong lithium absorption and NUV excesses of 2M1200- 
3405 and 2M1202-3328 confirm their youth and membership in 
TWA. Following the nomenclature for TWA members, we desig¬ 
nate 2M1200-3405 as TWA 35 and 2M1202-3328 as TWA 36. 
Neither star displays excess emission above photospheric levels in 
3-22 (im AllIV/SL photometry (Wright et al.|2010| >. The CMD po¬ 
sition of 2M105 8-2346 is consistent with a j3 Pic-like age or the 
upper envelope of Tuc-Hor members, but still marginally consis¬ 
tent with TWA. Until a better velocity measurement is available we 
refrain from assigning the star to either TWA or Octans. 


3.1.1 Possible binarity of2MASS J12002750-3405371 

The mean width of the cross correlation function (CCF) for 
2M1200-3405 (FWHM = 3.4 px) is much broader than the other 
candidates and standards (FWHM < 2.5 px; see fig. 4 of |Murphy| 
|& Lawson|20 15 1. This broadening is also visible in the raw spec¬ 
trum (Fig. [2} and Ha vio velocity width (Table[IJ, and is usually 


indicative of fast rotation (e.g. |Soderblom, Pendleton & Pallavicini] 
11989f >, unresolved spectroscopic binarity or a combination of the 
two. Using 2M1202-3328 as a narrow-lined template we can re¬ 
produce the level of broadening observed in 2M1200-3405 with 
either a vsin; or SB2 velocity shift of ~50 kms -1 , close to the 
cAA/A sa 45 km s~ 1 WiFeS velocity resolution. No broadening is 
evident at velocities <20 kms -1 (cAA/2A). Such a high implied 
rotation speed is unusual for single TWA stars (Jayawardhana et af| 
2006) and so may instead arise in the blending of lines from a close 
companion. However, in the absence of higher resolution observa¬ 
tions we cannot yet confirm this. If 2M1200-3405 is a binary then 
the unresolved CCF peak we measured should still trace the sys¬ 
temic velocity, assuming the binary is close to equal mass. 


3.1.2 2MASS J10284580-2830374 (=TWA 34) 


Another candidate, 2M1028-2830, is the recently-proposed TWA 
member and disc host TWA 34 (Schneider et al.||2012b| >. We 
confi rm the M4.9 spectral type and EW(Ha) = —9.6 A reported 
by Schneider et al. and measure EW(Li) = 600 mA and RV = 
20 ± 2 km s~ 1 from our higher-resolution WiFeS spectrum. From 
the modest strength and velocity width of its Ha line (viq = 
162 kms -1 ), TWA 34 does not appear to be accreting from its 
disc. However, we also detect strong forbidden [OI] A 6300 (EW = 
—3.0 A) and A6363 (—0.6 A) emission, indicating the presence of 


low-density gas usually attributed to a wind or outflow (e.g. Ap- 
penzellerjCMundt[ 1989[ |Rigliaco et al.|[2Q]~3] |Zuckerman, Vi- 
& Rodriguez (2014). Due to uncertainties in subtracting the 


strong auroral sky line, we do not unambiguously detect [Ol] 
A 5577 emission from any candidate. The SPM4 proper motion of 
TWA 34 yields a BANYAN II membership probability of 95 per 
cent, with an implied distance and radial velocity of 47 =t 6 pc 
and 13.5 ± 1.3 kms -1 . Rodriguez et al. (2015, in preparation) 
have obtained ALMA observations of TWA 34 and from a Kep- 
lerian fit to the CO velocity map derived a barycentric velocity of 
13.3 ±0.1 kms~*, in excellent agreement with the BANYAN pre¬ 
diction and other TWA members (Fig. [3}. This suggests that our 
WiFeS radial velocity may be erroneous, although we cannot yet 
rule out spectroscopic binarity. 

At 47 pc 2M1028-2830 is significantly under-luminous com¬ 
pared to other TWA members in Fig. [5] even considering reason¬ 
able errors on its photographic V magnitude. A distance of ~80 pc 
is required to bring it onto the TWA sequence of Riedel et al.| 
(2011) . [Schneider et al.|p0 12b) found congruent kinematic and 
photometric distances of ~50 pc when considering Mx and an em¬ 
pirical TWA isochrone, but the star is also under-luminous in that 
band. Its optical and near-infrared spectral energy distribution is 
well fit by a 3000 K model atmosphere (Bayo et al.|2008), consis- 
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Figure 6. Space motions of Sco-Cen members from |de Zeeuw et al.jj 1999] (red points) versus Galactic longitude (top) and latitude (bottom). Quadratic fits 
are shown to guide the eye. The mean velocities of the US, UCL and LCC subgroups (from left to right and top to bottom) are given by the blue triangles, as 
derived by|Chen et al. 1(271X1}. The space motions of 2M1239-5702 and 2M1422-3623 from Tableware plotted as large black points. 


tent with its V — K s colour and spectral type. A possible explanation 
for the underluminosity is that the photospheric light is intercepted 
by a flared disc seen at high inclination (similar to £ Cha 11: |Fangj 
|et al.|2013) . |Liu et al.|(2015) obtained Herschel/PACS photometry 
of 2M1028-2830/TWA 34 and modelled its disc. Although the best 
fit inclination angle of i = 60^j 5 deg is not strongly constrained by 
the data, their models rule out high inclinations (i > 80 deg). 


3.2 Two new accreting M-dwarfs in Sco-Cen 


The only two other stars in our sample with A11W75£ infrared ex¬ 
cesses, 2M1239-5702 and 2M1422-3623, lie within the classical 
boundaries of the LCC and and UCL subgroups of Sco-Cen, re¬ 
spectively (Fig. [TJ. Neither star is a member of the seven nearby 
moving groups tested by BANYAN (TWA, jS Pic, Argus, Tuc-Hor, 
Carina, Columba, AB Dor). The SPM4 proper motions and errors 
of 2M1239-5702 and 2M1422-3623 correspond to kinematic dis¬ 
tances of 119^[j pc and 142+24 pc, assuming the mean LCC and 
UCL space motions of |Chen et al.|(20TT) . These are in excellent 
agreement with the 120 and 140 pc mean distances found by [de| 
[Zeeuw et al.| ( |1999[ l for early and solar-type LCC and UCL mem¬ 
bers with Hipparcos parallaxes (also see Fig. |A 1 1 ). The CMD po¬ 
sitions of 2M1239-5702 and 2M1422-3623 in Fig. [5] suggest an 
age of approximately 10 Myr, similar to TWA and consistent with 
recent age ranges for UCL and LCC (Song et al. 2012, Peca ut et al.| 
|2012| >. The age of Sco-Cen and possible age spreads in the associ¬ 
ation are discussed in greater detail in Section |4~2| 

At these distances 2M1239-5702 and 2M1422-3623 are 
3 kms -1 and 6 kms~’, respectively, from the 


Chen et al. 


space motions (Fig. [ 3 }, compared to uncertainties of 2-4 km s -1 in 
each velocity component. However, unlike smaller, more coherent 
moving groups, OB associations are expected to show spatial and 
kinematic substructure. To examine the trend and spread of veloc¬ 
ities across Sco-Cen, we plot in Fig. [6] over 270 high-probability 
members from |de Zeeuw et al.| ( |1999p with radial velocities and 
updated Hipparcos astrometry. This is the largest available compi¬ 
lation of Sco-Cen members with six-dimensional phase space in¬ 
formation and a full description of its construction is given in Ap- 


pendix|A| The space motions of 2M1239-5702 and 2M1422-3623 
are perfectly consistent with the surrounding higher-mass members 
in this diagram. As expected, the mean subgroup velocities closely 
follow the longitude and latitude trends, albeit with significant scat¬ 
ter in individual velocities and minimal differentiation into the clas¬ 
sical subgroups. Much of this scatter likely results from observa¬ 
tional errors and single-epoch velocities of unknown spectroscopic 
binaries. For reference, |de Bruijne| ( |1999) found a one dimensional 
internal velocity dispersion of <1-1.5 kms~* from Hipparcos as¬ 
trometry. Figs. [6] and [AT] appear to confirm the suggestion of |Riz-| 
|zuto, Ireland & Robertson| $201 1j that Sco-Cen can be well mod¬ 
elled as continua in spatial position and velocity rather than three 
discrete subgroups. The addition of many more low mass members 
with radial velocities and Gaia parallaxes will better resolve the 
velocity and spatial substructure of Sco-Cen in coming years. 


The few Myr-old Lupus star forming region ( |Comer6n|2008| > 
lies within UCL and its members have similar distances and kine¬ 
matics to Sco-Cen. It is therefore plausible that 2M1422-3623 may 
be a young, outlying member of Lupus instead of the older Sco-Cen 
population. However, we do not believe this is likely. 2M1422- 
3623 is some 15-20 deg from the dark clouds which define the star 
forming region, and a recent kinematic and membership reanalysis 
of stars in the area by |Galli et al.| l |2013j l did not find any legitimate 
‘off-cloud’ Lupus members within 9 deg of 2M1422-3623. The 
results of |Galli et al.|(20l3| strongly support the hypothesis that 
the majority of the ~150 late-type, older pre-main sequence stars 
found around the Lupus clouds in X-ray surveys (e.g. |Krautter et al.] 
| 1997|[Wichmann et al. 1 997} are in fact members of UCL/Sco-Cen 
l |Mamajek et al.|2002| | Aeibisch & Mamajek|2008| l. 2M1239-5702 
is situated close to the centre of LCC, only 1 deg from the fore¬ 
ground M-giant and Southern Cross member y Cru. It lies a few 
degrees north of the so-called ‘Park & Finley’ stars (|Park & Finley| 
|1996[ |Preibisch & Mamajek|2008| >, a group of four young K and 
M-type LCC members discovered in the vicinity of the B0.5 mem¬ 
ber j8 Cru l |Feigelson & Lawson|1997||Alcala et al.|2002fr . None of 
these stars appear to be accreting or host circumstellar discs. 
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Figure 7. Spectral energy distributions (SEDs) of 2M1422-3623 (left) and 2M1239—5702 (right). Photometric measurements from DENIS (/), 2MASS and 
WISE are plotted in red. The total SED (blue line) is approximated by a solar metallicity BT-Settl model (grey line; |Allard et al.|2012fr fitted to the DENIS 
and 2MASS data, and two blackbodies (green and orange lines). The characteristic emitting area of each component is given by the factor T = Lir/L*. The 
interquartile range of K5-M2 Class I Taurus sources from[D’Alessio et al.|{1999) is shown by the shaded region, normalised at 1.6 pm. 


3.2.1 Disc and accretion properties 


The spectral energy distributions of 2M1422-3623 and 2M1239- 
5702 are compiled in Fig. [ 7 ] Both stars exhibit strong excesses in 
the W3 (12 pm) and W4 (22 pm) WISE bands. These are con¬ 
sistent with evolved, ‘homologously depleted’ discs which are be¬ 
coming optically thin without forming an inner hole or gap, and 
show a strong flux decrement at all infrared wavelengths relative to 
optically-thick ‘full’ discs (Espaillat et al. 2012, Luhman & Mama- 
|jek|2QI2} . The observed (K s — W4) excesses of 2M1422-3623 and 
2M1239-5702 are similar to discs around K5-M5 stars in Upper 
Scorpius (age ~10 Myr), which already show evidence of signif¬ 
icant evolution compared to Taurus ( |Luhman & Mamajek|2012} . 
Although not necessarily physical, double-temperature blackbody 
models with a cool ~230 K outer component and a warmer 800- 
900 K inner component provide a good fit to the WISE photometry 
in both cases. Assuming L* = 0.03 Lq (Fig.JTjl and large (black- 
body) grains, these temperatures correspond to disc radii of 0.25 au 
and ~0.02 au, respectively (Backman & Paresce 1993). Further 
infrared and millimetre observations (e.g. with VLT/VISIR and 
ALMA) will help elucidate the structure and chemistry of the discs 
around 2M1422-3623 and 2M1239-5702, for instance by resolv¬ 
ing the 10 pm silicate feature if it exists, or detecting cold dust and 
CO gas emission. 

Accretion of gas from the inner disc onto the star is usually 
accompanied by enhanced Balmer and other line emission (|Muze-| 
|rolle, Hartmann & Calvet|1998} . Fla velocity profiles for 2M1239- 
5702 and 2M1422-3623 are given in Fig. [8] The broad and strong 
emission at all epochs (also see Table|3j confirms both stars are 
accreting, with EWs greater than the 15-18 A threshold for M3- 
M6 stars suggested by Fang et al. (2009jl and vp) velocity widths 
in exces s of the 270 kms 1 accretion criterion of lWhite & Basril 
(2003^ Both stars presented weaker Fla emission in the second- 
epoch observations, falling just below the 270 km s _1 threshold but 
still above the EW limit. Subsequent observations of 2M1422-3623 
show a recovery in both EW and vio- Applying the vp)(Fla)-M 
relation of Natta et al. (2004| > to the velocities in Table [3] yields 
mass accretion rates of ~10~ 9,5 Mq yr~* at peak and a factor 
of ten lower at minima, comparable to those in older groups like 
TWA and tj Cha i jMuzerolle et al.|[2000| [Lawson, Lyo & Muze-| 
|rolle||2004j>. Similar behaviour is seen in the ~6 Myr-old rj Cha 


4 Jayawardhana, Mohanty & Basri 1 2 003} suggest a lower accretion thresh¬ 
old of 200 km s for M5-M8 objects. Whatever the adopted criterion, it is 
clear from the large EWs and the presence of other indicators that 2M1239- 
5702 and 2M1422—3623 are accreting from gas-rich inner discs. 



Figure 8 . 2.3-m/WiFeS RV-corrected Ha velocity profiles of the accretors 
2M1422-3623 (top) and 2M1239-5702 (bottom). For comparison, each 
panel also shows the profile of the new TWA member and non-accretor 
2M1202-3328 (shaded region; EW = -9.5 A, v 10 = 169 km s^ 1 ). 


members RECX 5 and 2M0820-8003, which are borderline accre¬ 
tors but have undergone sporadic episodes of strong accretion in 
multi-epoch observations (Lawson et al. 2004, Murphy et al.|20lT} . 
Unlike RECX 5, which hosts a transitional disc with an inner hole 
and gap (Bouwman et al.|2010} , 2M0820-8003,2M1239-5702 and 
2M1422-3623 show excesses at <5 pm, indicating the presence of 
dust in the inner disc. Their variable accretion is presumably driven 
by ‘clumpy’ residual gas-rich material in this region. 

In addition to Ha we also observed Nal D and forbid¬ 
den [Ol] A6300/6363 emission in 2M1422-3623 and strong He I 
A5876/6678 emission in both stars (Fig. [ 7 ] and Table[3}. The 
strength of the Na I and He I emission was generally correlated with 
Ha, as expected from their common origin in magnetospheric infall 
regions i jMuzerolle et al.|l998| , whereas the outflow-driven forbid¬ 
den emission remained approximately constant. Continuum excess 
is another indicator of accretion and activity in T Tauri stars, partic¬ 
ularly in the UV ( [Bertout, Basri & Bouvier| 1988 j i. 2M1239-5702 
and 2M1422-3623 both show enhanced GALEX near-UV emission 
compared to the other lithium-rich stars and most stars in the Torres] 
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8 Murphy, Lawson & Bento 


Table^ ; J3missionUne i ec^iivalentwidthsJn2M1239^5702^md2M1422^3623 ; _Thevelocit^w^ 

Estimated uncertainties are ±10 km s 1 for i ! jq, +0.1 A for all lines other than Ha and ±2 A for the broader Ha line. 


Star 

UT Date 

vio(Ha) 

(kms -1 ) 

Ha 

(A) 

He I A5876 
(A) 

He IA6678 
(A) 

[Ol] A 6300 
(A) 

[Oi] A 6363 
(A) 

Nal A5896 
(A) 

Nal A5890 
(A) 

2M1239-5702 

2015 March 9 

331 

-63 

-5.3 

-1.2 



-0.2 

-0.2 


2015 March 30 

238 

-27 

-2.4 

-0.4 

-0.4 


-0.6 

-0.9 

2M1422-3623 

2015 March 8 

335 

-91 

-5.7 

-1.3 

-2.5 

-0.5 

-3.2 

-1.9 


2015 March 31 

236 

-33 

-3.4 

-0.6 

-2.8 

-0.8 

-2.0 

-1.6 


2015 April 1 

303 

-54 

-1.7 

-0.3 

-3.7 

-1.0 

-1.8 

-1.0 


2015 April 6 

341 

-68 

-3.1 

-0.9 

-2.4 

-0.6 

-2.4 

-0.8 


jet al.|p008) sample of young moving group members. Comparing 
Li I A6708 line strengths we deduce that only the 2015 March 8 
spectrum of 2M1422-3623 was noticeably veiled (when accretion 
was strongest). On that night we measure an EW of 550 mA, com¬ 
pared to values of 600-630 mA at the other three epochs. 


3.3 Other lithium-rich stars 


2M1005-2131 and 2M1049-2509 both have smaller proper mo¬ 
tions, inconsistent with neighbouring TWA members and indicative 
of larger distances. We find satisfactory matches to the TWA space 
motion at 75 and 95 pc, respectively (Fig. [3}. However, at these dis¬ 
tances the stars are very far from other TWA members in space (e.g. 
jWeinberger et al.||2013| l and still sit below the mean CMD trend 
in Fig. [5] Adopting the UCAC4 proper motion for 2M1005-2131 
(—55.3±7.1, — 9. 7±5.3 mas yr _1 ) increases its BANYAN member¬ 
ship probability from 13 to 50 per cent (including RV), but at the 
implied distance of 42 pc the star lies near the Pleiades sequence 
and its predicted proper motion is much larger than observed. 

2M1049-2509 was observed with VLT/UVES on 2014 
May 11 (programme ID 093.C-0133A) as part of GALNYSS. Ex¬ 
amining the pipeline-reduced archival spectrum we find EW(Li) = 
680 ±30 mA, EW(Ha) = -22 ± 3 A, and -300 mA Call in¬ 
frared triplet emission. The radial velocity of the Nal A8200 ab¬ 
sorption doublet (17± 1 kms -1 ) agrees with our lower-resolution 
WiFeS measurement. The star was also serendipitously observed 
by both the Chandra and XMM-Newton X-ray satellite s. Its 0.2- 
12 keV XMM luminosity at 95 pc (10 28 9 erg s~ 1 ; Rosen etal.| 2015) 


and log Lx /±bol ~ —3.0 are both consistent with an age less than 
the Pleiades. Both 2M1005-2131 and 2M1049-2509 are clearly 
young, but their exact membership awaits further observations, par¬ 
ticularly parallaxes. 


4 DISCUSSION 

4.1 Other accretors in Sco-Cen 


Outside of Upper Scorpius, 2M1239-5702 and 2M1422-3623 join 
only a handful of known accretors in the older UCL and LCC sub¬ 
groups. jS^MTn^eraL] l|2(TITjl suggested the M3.5 LCC member 
2MASS J11311482^-826279 was accreting based on its broad Ha 
line (vio = 233 kms -1 ). However, given its weak EW (—7.3 A) 
and that the star has no excess emission visible in WISE, we sug¬ 
gest the broad Ha line |Shkolnik et al.[ observed was the result 
of chromospheric activity. Flares can give rise to broad, multi- 
component Ha emission like that seen in 2M1131^4-826 and are 
usually characterised by lower EWs than observed in accretors 
JMont£»et^aL 1998| [Jayawardhana et al,|[2006| jMurphy et al.j 
[2011 1 . Rodriguez et al. proposed another M3.5 LCC accre- 

tor, 2MASS J13373839-4736297, based on a strong (-13.7 A) and 
asymmetric Ha line. jSchneider et al.j ( [2012bj l confirmed the star 
hosts a disc with excesses in all four WISE bands. Accretion thus 


seems a likely explanation for the observed line emission. |Zucker-| 
man]( |2015) report 2M1337-4736 is a 10 arcsec wide binary of near 
equal brightness. The secondary does not appear to possess a disc. 

The three confirmed M-type accretors in Sco-Cen (2M1337- 
4736A, 2M1239-5702 and 2M1422-3623) are joined by the Her- 
big Ae star HD 139614 and two accreting F stars; AK Sco and 
HD 135344 fPreibisch & Mamajek|j2008j| l. These authors noted 
HD 139614 and 135344 lie near the western edge of the Lu¬ 
pus clouds and so may be younger than the median age of UCL 
(16 Myr; jPecaut et al.||2012] >. jMamajek et al.j ( [2002) proposed 
the K1 T Tauri star MP Mus (=PDS 66) as a southern member 
of LCC. It was the only star which showed signs of accretion 
(EW(Ha) = -39 A) in their survey of 110 G and K-type UCL 
and LCC stars. However, jMurphy et al.j ( |2013[ > have disputed the 
membership of MP Mus in LCC, claiming a better kinematic and 
distance match to the 3-5 Myr e Chamaeleontis association which 
overlaps the south of LCC. 

In light of the increasing number of accretors in Sco-Cen, 
membership of the inconclusive TWA member and confirmed ac- 


cretor 2MASS 112071089-3230537 (TWA 31; [Shkolnik et ak| 
|201 1| jSchneider et al. 2012a) deserves reevaluation. TWA 31 
was rejected by Ducourant et aL] ( j20l4) on the basis of its poor 
USNO-Bl.O proper motion. BANYAN gives a membership proba¬ 
bility of only 28 per cent at 59 ± 7 pc. Its SPM4 proper motion, 
(—52 ± 8, —18 ± 8) masyr -1 , is a better match to TWA, with a 
membership probability of 98 per cent at 61 ± 7 pc. However, with 
Vspm4 = 18.45 (photographic) TWA 31 would be under-luminous 
at this distance. It is not found in UCAC catalogues and its PPMXL 
proper motion agrees with USNO-B1.0. At its proposed photomet¬ 
ric distance of 110 pc i jShkolnik et al.|20TT) , TWA 31 has an SPM4 
space motion somewhat similar to other TWA members and a rea¬ 
sonable CMD position, but is much further away than the majority 
of TWA stars. Only TWA 15AB is more distant (d ~ 115 pc), but 
it is far enough south to conceivably be an LCC member (but see 
discussion in jDucourant et aL]|2014j ). Until better astrometry are 
available the final status of TWA 31 remains unclear. 


4.2 Implications for the age of Sco-Cen and disc lifetimes 


The fraction of accreting stars in young groups and star forming 
regions is observed to fall rapidly with age, with e-folding time- 
scales of 2-3 Myr i Mamajek||2009| jFedele et al.|j2010] >. For in¬ 
stance, jFrasca et al, 1 2015) recently found accretor fractions of 
35-40 per cent for the ~2 Myr-old Cha I population, falling to 
only 2^4 per cent in the older 010 Myr) y Velorum cluster. Sim¬ 
ilar behaviour is observed in young moving groups, where the ac¬ 
cretor fraction in the <10 Myr q and e Cha groups is ~30 per 
cent i jMurphy et al. ||2013) but by the age of Upper Scorpius (10 
Myr) has fallen to ~10 per cent ( jLodieu et al.|20lT|) and in UCL 
and LCC (16-17 Myr) is at most a few per cent ( |Mamajek et al.j 
|2002) . Complicating matters are recent recalibrations of pre-main 
sequence cluster and association ages by|Bell et al.j (20 1^3), which 
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Figure 9. Li I A6708 EWs of 2M1239-5702 (star) and 2M1422-3623 (tri¬ 
angle) compared to new UCL and LCC members from Song et a~T~|^20 ] 2} 
Members of TWA ^Schneider et al.|2012a) and j8 Pic jda Silva et al'|2009| 
|Binks & Jeffries 2014 1 are plotted for comparison. The —20 Myr-old j8 Pic 
LDB is also given (dashed line). Sco-Cen members have lithium depletion 
levels between TWA and Pic, corresponding to an age of 10-20 Myr. 
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have increased the ages of some benchmark groups by up to a factor 
of two and therefore increased average circumstellar disc lifetimes. 
Nevertheless, if 2M1239-5702 and 2M1422-3623 are at the me¬ 
dian age of LCC and UCL then they are incredibly rare objects. As¬ 
suming exponential decay is an appropriate functional form for the 
di ssip ation of inner discs at older ages (see discussion in|Mamajek| 


2009), for a putative population of 10 4 Sco-Cen members we may 


(very) naively expect only 5-10 stars to show signs of accretion at 
16-17 Myr given an e-folding time-scale of z = 2.3 Myr ( |Fedele[ 
|et al.|20K)} . Even with accretor fractions of 1-2 per cent (consis¬ 
tent with t ~ 8-10 Myr) there should exist no more than a few 
hundred examples across the association. The serendipitous discov¬ 
ery of two such stars in a kinematically-selected survey not look¬ 
ing at Sco-Cen or targeting accretors and infrared excess sources 
thus seems very surprising, and demands a more nuanced expla¬ 
nation than simple exponential decay at fixed age and decay time. 
In the following paragraphs we consider two possible solutions to 
this problem - significant age spreads within Sco-Cen and stellar 
mass-dependent disc dissipation. 

After considering the effects of observational uncertainties, 
|Mamajek et al.|p002| ) concluded that the intrinsic 1 a age spreads 
of G-type stars in UCL and LCC are 3 Myr and 2 Myr, respec¬ 
tively. This implies that almost all star formation ceased in the sub¬ 
groups approximately 5-10 Myr ago. |Pecaut et al.| ) |2012j ) found 
similar median ages for UCL (16 Myr) and LCC (17 Myr) F-type 
members, but with larger intrinsic age spreads of ±4-7 Myr and 
±0-8 Myr, depending on the choice of evolutionary models. In the 
maximal cases these give 2a limits on the cessation of star forma¬ 
tion which are comparable to the ages of the subgroups themselves. 
Such age spreads are not unexpected given that star formation in 
UCL and LCC probably unfolded in a series of small embedded 
clusters and filaments in a large, dynamically unbound complex 
l |Preibisch & Mamajek[2008$ . Given the large age spreads observed 
in higher-mass Sco-Cen members, it is therefore plausible that 
2M1239-5702 and 2M1422-3623 may be as young as 5-10 Myr 
and among the last generation of stars to form across the region. 
In addition to ongoing accretion from dusty inner discs, young 
ages for 2M1239-5702 and 2M1422-3623 are also supported by 


their elevated CMD positions and undepleted photospheric lithium 
(Fig.0. 

A second possibility is that the time-scale for circumstellar 
disc dispersal is stellar mass dependent, with the processes respon¬ 
sible for clearing optically-thick primordial discs being less ef¬ 
ficient around lower-mass stars (e.g. |Kennedy & Kenyon||2009| l. 
Such a trend has been observed in several young groups, including 
Upper Scorpius (Carp enter et al.|2 006; Luhman & Mamajek 2012| >. 
These studies indicate that a significant fraction (~25 per cent) of 
inner discs around low-mass stars can survive for at least 10 Myr. 
There are hints such a phenomenon may extend to the rest of Sco- 
Cen. Modulo the association mass function, three M-type accretors 
have now been identified in UCL and LCC, the same number of ear¬ 
lier type members showing signs of accretion. Larger numbers of 
K and M-type members are necessary to test this hypothesis. Both 
2M1239-5702 and 2M1422-3623 have strong near-UV excesses 
and were included in the spectroscopic sample for this reason. The 
inclusion of GALEX photometry in future Sco-Cen surveys should 
be a useful way to select additional active and accreting objects 
(e.g. |Rodriguez et al.|201l||2013||Shkolnik et al.|201l) . 

Finally, we comment on the recent suggestion of |Song et al.| 
( |2012j ) that UCL and LCC are as young as ~10 Myr. This age esti¬ 
mate was based on ~90 GKM stars which showed levels of lithium 
depletion bracketed between TWA (8-10 Myr) and the /3 Pictoris 
Association (Fig. [9}, which had previously been considered to have 
an age of ~12 Myr. Given recent upward revisions of the age of 
)8 Pic to 20-25 Myr |Binks & Jeffries||2014| |Malo et al.||2014|), 
the |Song et al,| ( |2012 i results now imply ages between 10 and 
20 Myr, consistent with the 16-17 Myr median subgroup ages and 
age spreads presented by |Pecaut et al.|p() 1 2) and discussed in this 
section. Again, many more M-type Sco-Cen members are needed 
to confirm the lithium depletion age of the association and examine 
trends with stellar parameters such as rotation and binarity. 
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APPENDIX A: SCO-CEN MEMBERS WITH HIPPARCOS 
ASTROMETRY AND RADIAL VELOCITIES 

|Rizzuto et al.| ( j2011[) presented a list of 436 Sco-Cen members 
from the new Hipparcos reduction (van Le euwen 20071, includ¬ 
ing 382 stars with radial velocities in the CRVAD-2 compilation 
( [Kharchenko et al.|2007ft , which they used to construct linear mod¬ 
els of the Sco-Cen space motion with Galactic longitude. During 
the present work we discovered that most of the radial velocities 
in that sample are non-spectroscopic ‘astrometric’ velocities from 
|Madsen, Dravins & Lindegren| ( |2002] l which are inappropriate for 
calculating space motions or memberships. These were apparently 
sourced from the Pulkovo Compilation of Radial Velocities (PCRV; 
|Gontcharov||2006| , however they are not in the published PCRV 
catalogue. They were part of a preliminary list used to construct 
CRVAD-2 and later removed prior to publication of the PCRV 
(G. A. Gontcharov, personal communication). 

To construct Fig. [6] we therefore assembled our own com¬ 
pilation of purely spectroscopic velocities from several large 
catalogues in the literature. The sample was constructed as follows. 
We first crossmatched the 521 members proposed by |de Zeeuwj 
|et al.| ( |1999| ) (dZ99) against |van Leeuwen| [2007[ l and the list 
of Sco-Cen stars with Magellan/MIKE velocities from |Chen| 
[et al.| p01 1} . This returned 181 matches. After removing known 
spectroscopic binaries and non-members 140 stars were retained. 
We chose as a secondary source the PCRV, which contains mean 
absolute velocities for 35,495 Hipparcos stars from over 200 
publications, including major surveys. This returned 134 matches 
against dZ99, and we included the 110 stars not already observed 
by jChen et al.| Finally, we crossmatched dZ99 against CRVAD-2, 
after removing all objects sourced from the pre-publication PCRV 
list (ICRV=2). This gave 57 primaries not already in the PCRV 
or |Chen et al.| samples. Several stars were found in more than 
one catalogue, allowing us to eliminate the spectroscopic binaries 
HIP 82747 and HIP 73666. The result of these crossmatches 
was a list of 305 Sco-Cen members with high quality astrometry 
and radial velocities (median precision 1.5 kms~'). The final 
compilation is listed in Table | A1 1 and plotted in Fig |A 1 [ Note that 
Fig. E includes additional cuts to remove 31 stars with unusually 
large velocities (—20 kms -1 > RV > 25 kms -1 ) or distances 
0 d > 200 pc) which are less likely to be members. 


SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online 
version of this article: 

Table |Al| New Hipparcos astrometry, literature radial velocities 
and re-calculated UVW space motions and XYZ positions of 305 
Sco-Cen members proposed by de Zeeuw et al.|(T999[ l. 

Please note: Oxford University Press are not responsible for the 
content or functionality of any supporting materials supplied by 
the authors. Any queries (other than missing material) should be 
directed to the corresponding author for the article. 

This paper has been typeset from a TpX/I6Tj:X file prepared by the author. 
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Table Al. Radial velocities and space motions of Sco-Cen members proposed by |de Zeeuw et al.]{1999) from the Hipparcos catalogue. The final column lists 
the source of the radial velocity; either [Qien et al.|{2011) , PCRV (Gontcharov|2QQ6| or CRVAD-2 (Kharchenko et al.|2007) . The full table with extra columns 
is available in machine-readable format as Supporting Information with the online version of the article and on the CDS VizieR service. 
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Figure Al. Position, radial velocity and parallax distributions of the stars in Table |Al| The source of the RV is indicated: |Chen et al. 1(2011) (red diamonds), 
PCRV (blue squares) or CRVAD2 (green dots). Crosses show the 31 potentially unreliable stars not included in Fig. [6] These have d > 200 pc or —20 km s -1 > 
RV > 25 kms -1 (dashed lines). 
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